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Abstract 


We use the epsilon expansion to explore a new universality class of second order 
quantum phase transitions associated with a four-dimensional Yukawa field theory 
coupled to a traceless Hermitean matrix scalar field. We argue that this class includes 
four-fermi models in 2 < D < 4 dimensions with SU(Np) x U(N ) symmetry and a 
U(N ) scalar, SU(Np) iso-vector 4-fermi coupling. The epsilon expansion indicates that 
there is a second order phase transition for N > N*(Np), where N*(Np) — .27 Np if 
Np —■* oo. 
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There are several scenarios in field theory where quantum critical phenomena play a role. 
They are intimately associated with renormalization theory and are important in modifica¬ 
tions of the standard model using technicolor U or models with composite Higgs fields @- 
[6J. Quantum critical phenomena are associated with phase transitions which are achieved, 
typically at zero temperature, by tuning mechanical parameters such as particle masses or 
coupling constants to critical values. Classical critical phenomena, on the other hand, occur 
at non-zero temperature and are associated with phase transitions where temperature or 
some other thermodynamic variable is the quantity which must be tuned. 

An important concept in the study of second order phase transitions is universality. 
A phase transition is characterized by the number of degrees of freedom which become 
massless at the critical point, as well as the symmetry of the system and the symmetry 
breaking pattern associated with the transition. In classical critical phenomena, at finite 
temperature, fermions decouple from infrared behavior because of their half-odd-integer 
Matsubara frequency. Thus, all universality classes are represented by Bose field theories. 
In quantum critical behavior, fermions can be massless at a phase transition and therefore 
they must be included in the effective low energy field theories which represent universality 
classes. 

A useful tool for the study of phase transitions is the Wilson-Fisher epsilon expansion 
where the beta function is computed in a 4 — e dimensional effective field theory and second 
order phase transitions are characterized by the infrared stable fixed points of the renor¬ 
malization group flow. If the beta function exhibits such fixed points and if the coupling 
constants have initial conditions within their domain of attraction, there is a second order 
transition. If no infrared stable fixed points exist then if there is a phase transition it must 
be of first order, driven by fluctuations 0-[0. 

The epsilon expansion has been used to show that the 0(N ) vector model in 4 — e di¬ 


mensions has a second order phase transition with computable critical exponents |L0|. It 
represents a universality class which contains the O(N) nonlinear sigma model in dimen¬ 


sions between 2 and 4. The sigma model is renormalizable only in the 1/N expansion [11 


and represents the continuum limit of generalized magnetic systems. It has also been ar¬ 
gued |12[]-[[T4|] that the four dimensional Yukawa-0 4 theory in 4 — e dimensions represents 
the universality class of the Gross-Neveu model in dimensions between 2 and 4. The latter 
is also renormalizable only in the 1/N expansion []T3], [RJ and the critical exponents have 
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been computed to high order |T7|-||T9|]. The 0(N ) model has one relevant coupling constant 
in the 4-dimensional theory, so the fixed point of order e is automatically infrared stable 
and the phase transition is of second order. When there is more than one relevant interac¬ 
tion in the 4 — e dimensional theory, and therefore a multi-dimensional coupling constant 
flow, the generic situation is that fixed points of the flow are not infrared stable and the 
phase transitions are first order. The classic example of this is the massless scalar electro¬ 
dynamics studied by Coleman and Weinberg [Q. There, in the two-dimensional plane of 
electromagnetic and Higgs self-coupling there are no infrared stable fixed points. Another 
example occurs in the complex matrix scalar field theory with SU(N F ) x SU(Np) symmetry 
which represents the universality class of the finite temperature chiral transition in QCD 
|20| . There are two coupling constants for the renormalizable interactions tr(ibFM) 2 and 
(trM^M) 2 and when the matrices are larger than 2x2 the phase transition is first order. 
When the matrices are 2x2 the model is equivalent to a vector theory with one coupling 
constant and the phase transition is second order. Similar reasoning has been used to ar¬ 
gue that a 4 — e-dimensional traceless Hermitean matrix scalar field theory with SU(Np) 
symmetry, which also has two relevant couplings tr</> 4 and (tr</> 2 ) 2 , has a second order phase 
transition only when Np = 2 and has a fluctuation induced first order phase transition when 
Np > 2 |2T| . 

In this Letter, we shall show that when a Yukawa coupling to fermions is introduced 
to the Hermitean matrix theory, the existence of a nontrivial Yukawa fixed point tends to 
make the fixed points of the matrix theory infrared stable. For a given Np, we always find a 
second order phase transition if the number of fermions in the fundamental representation is 
large enough, asymptotically N > .27 Np for large N, Np, and first order otherwise. We shall 
also argue that the Yukawa-matrix model represents the universality class of a second order 
quantum phase transition in a matrix-coupled four-fermi theory in 2 < D < 4 dimensions. 
We use this equivalence to argue that, for large enough N, the epsilon expansion is valid for 
epsilon of order 1. 

The Euclidean action for the 4 — e-dimensional Yukawa model is 


S = 


d“-‘x I « I V? - V + ) </>? 


1 -=>,=> 8n 2 /j, e 

-tr F V</> ■ V</> + — 


M (t rW,v + 


(i) 
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where a, b — 1,... N, a, (3 — 1,..., Np, 0 Qj g is a traceless Hermitean Np x Np scalar matrix 
field. -0“ is a complex Dirac spinor. Here we have included all vertices which are compatible 
with the unitary symmetry, 0 —> UcftU^, 0 —> 170, 0 —> ipU' with 17 G SU(Np) and which 
are renormalizable in 4 dimensions. We have scaled the coupling constants so that in the 
large Np limit planar graphs dominate and fermion loops are suppressed and in the large N 
limit the bubble graphs dominate. 

The beta functions of this model for the case y — 0 was computed to one-loop order in 
The beta function for the model with y ^ 0 is readily obtained to the same order: 
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■y 2 gi, 


9 
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8NN F 


y A + 


2 N, 


-y 2 92 , 


e , 2NNp + N 2 F -3 3 

P v = ~y + — 1gA r2 at —y 


T ' 16NjN 
Also at one-loop order, the scaling dimension of the scalar field is 


( 2 ) 

(3) 

(4) 


A b — 


8 N, 


-r + •■ ■, 


(5) 


and that of the fermion field is 


Af = 


3 — e 


1 N 2 F 


32 NpN 


y + • • • 


( 6 ) 


There will be a second order phase transition when the beta function has infrared stable 
fixed points. This occurs where A A) = 0 and where all eigenvalues of the matrix dfii/dgj 
are positive. The trivial fixed point y^ = 0 of the Yukawa coupling always coincides with a 
negative eigenvalue of the stability matrix. It could therefore be ultraviolet stable, but not 
infrared stable. Thus, to find infrared fixed points, we take the no n zero fixed point of the 
Yukawa coupling, 

RN'IN 

:e- (7) 


yl = 


8 N 2 N 


2 NN f + JVj. - 3' 

This coupling is small when e is small. 

In the large N limit, the beta functions have an infrared stable fixed point at Yukawa 
coupling (J7j) and scalar self-couplings 


9 1 = 


18 {N 2 + 3) 
N 2 


* _ 6 Np 

e + ... , g 2 - —e + ... 


( 8 ) 
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In this regime, the scalar self couplings are small and the one-loop perturbation theory is 
accurate. The scaling dimensions of the spinor and scalar fields can be computed to leading 
order in e and 1/A at the fixed points 


Ap = 


e + N l 


2 8 N Np 
A b = 1 + ... 


■e + ..., 


(9) 

( 10 ) 


These critical exponents depend on both A and Np and therefore each value of these pa¬ 
rameters specifies a different universality class. 

We have solved for the fixed points numerically for general values of A and Np. The 
asymptotic form of the fixed points in (0) and © continued to lower A are depicted for a few 
values of Np in Fig. 1. A critical value of A, N*(Np), can be defined for every Np, such that 
an infrared stable fixed point exists at A > N*(Np). If Np —> oo then N*(Np) ~ 0.27Np. 
This asymptotic limit is reached through values of IV*(2) = 34, A*(4) = 8, A*(8) = 4, 
A*(80) = 22, A*(800) = 216. 

An example of a dynamical system whose universality class is described by © is the 
four-fcrmi theory [^2[ 


S = 


d°x i rj ■ - p I rj'£ l3 r f) r 1 T«M 


A 


Aj.b 


A 


( 11 ) 


where 2 < D < 4, a, /?,... = 1,..., Np, a,b } ... = 1,..., A and ip is a 4-component Dirac 
spinor. T A are generators of SU(Np) in the fundamental representation, normalized so 
that trT A T B = 5 AB /2. This model has the same symmetries as (0): C, P, T, as well as 
discrete chiral symmetry and global SU(Np) x 17(A). [^3| The parity and discrete chiral 
transformation of the four component fermions can be defined in such a way that the mass 
operator ijnp transforms as a scalar under parity and changes sign under a discrete chiral 
transformation [^4J] . The phase transition that we shall study generates fermion mass through 
spontaneous breaking of the SU(Np) and discrete chiral symmetry. The order parameter is 
0a/3 = ( — A/A) < 1p“lp a p - -^5 af3 1p“lp“ >. 

Four-fermi models with isoscalar couplings are known to be renormalizablc in the large A 
I. pR |26| when the coupling is tuned to a critical value where a second order 


expansion 


15 


phase transition takes place. Existence of a phase transition has also been demonstrated in 


the conformal bootstrap approach in isoscalar models [27|. Their arguments can be extended 
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to the model (|lTD in a straightforward way. We introduce an Np x Np auxiliary held, 0, 

S = f d D x ($“{S aP 7 • V + 0 a/ 3 )0^ + ^trp0 2 ^ i ( 12 ) 

where tr^0 = 0. Here, the fermions can be integrated to obtain the scalar held theory 

S eff = -NTR In (7 • V + 0 ) + J d D x^tr F <j) 2 , (13) 

which can be studied perturbatively in the 1/N expansion. (Here TR means trace in function 
space and over indices.) The effective scalar propagator for the large N expansion, 

^a/3ys(p) — 4 > a(d(p')4 > 'ys(. p) ^ ^aS^/3 7 jy 

x 1 f 1 1 4r(l-D/2)\r(D/2)]*p D -* \- 1 

N\ A A c (Att) d / 2 T(D -1) y ■ 1 J 

A c is the critical coupling which is of order A 2 ~ D where A is the large momentum cutoff. The 
propagator is scale invariant when A is tuned to A c . This is the value of the coupling at which 
the phase transition takes place. When A has this critical value, the scaling dimensions of the 
scalar and spinor helds are defined by the behavior of the fermion and scalar propagators, 


Sapip) ~ N 2Af+1 d , (15) 

Aa^sip) ~ (SasSfa - j^8 a p8 lS ) b! 2As_D , (16) 


respectively. From this and computation of the spinor self-energy to leading order in 1/N 
and 0), we deduce 


D- 1 T(D — l)(Np — 1) 

— 2 8 N^NT(i - D/2)[r(D/2)] 2 r(D/2 - 1) + 


(17) 


A B = 1 + ... 


(18) 


If we set D = 4 — e in (ITT]) and expand to first order in e, the scaling dimension of the fermion 
in the four-fermi model is identical with that in the Yukawa-matrix model at the infrared 
stable fixed point in (||). Furthermore, the scaling dimension of the scalar held to leading 
order in (|l^) agrees with that in (|T 0 |). We conjecture that this equivalence holds in higher 
orders of perturbation theory and that the two models are in the same universality class. 

Symmetry breaking in the model m can also be studied in the 1/N expansion by 
computing the effective potential for 0 = (0) which to leading order in N is readily obtained 
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from (p~3| ) as 


J5S. (A = (i - £) 5*^ + 2 1^ tr 'C) D/2 ' < 19 > 

This potential clearly exhibits a second order phase transition at A —> A c- When A > A c, 
the eigenvalues of 0 obtain an expectation value with equal magnitude and undetermined 
sign. This symmetry breaking is consistent with the fact that tuj) = 0 only if N F is even 
and the symmetry breaking pattern is SU(N F ) —■> SU(N F / 2) 0 SU(N F /2). (If N F is odd, 
we can show for small values of N F that the symmetry breaking pattern with [N F — l)/2 
positive, (N f — l)/2 negative and one zero eigenvalue is preferred. In the following we shall 
avoid this complication by assuming that N F is even.) If we assume this symmetry breaking 
pattern, it is straightforward to compute the leading correction to (^) when D = 3 and 
A = A c- The 1-loop correction to the effective action is given by [^2| 


I4r(0) = —A+W|0| 3 


4 AT-2 


In 101/// + ...) NN F \<f>\ 


„ 2 N f -2 

1 V 


( 20 ) 


7T 2 N 

The logarithmic corrections which are typical of next to leading order in the effective potential 
(and which can in some cases lead to a Coleman-Weinberg instability ) exponentiate to give 
the scaling dimension of the potential (22)]. Thus we have established that if N is large 
enough, the model has a second order phase transition. For large N F the critical number of 
colors is N* = ,27N F . At N F = 2 N* = 34, at N F = 4 N* = 8, at N F = 200 N* = 55. The 
behavior of the critical number of colors as a function of the number of flavors is consistent 


with results obtained in equivalent four fcrmi interaction theories p2|. It is interesting to 
speculate that the upper critical value of N F ~ 3.7 N for existence of a chiral transition in 
the four-fermi theory is related to the upper critical N F ~ const. • iV co i or for existence of 
chiral symmetry breaking in 2+1 dimensional QCD |^D[. A correspondence of this kind was 
the intention of ref. |^TJ . However, if we naively identify U(N) with color, that theory has 
either different symmetry or different degrees of freedom: If the QCD is confining, U(N) is 
absent and if it is not confining, massless gluons should contribute to the critical behavior. 
Furthermore, in gauge theories at zero temperature, once N F and N are given, either the 
vacuum is symmetric or the symmetry is spontaneously broken. There is no mechanical 
parameter that drives the symmetry breakdown unlike at finite temperature. Suppose an 
additional interaction such as the four-fermi coupling that we have considered in this paper is 
included in the theory. Then this type of effective low energy field theory might be applicable 
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for the gauged Nambu-Jona Lasinio model and our analysis would apply close to the Nambu- 
Jona-Lasinio fixed point where the gauge interaction is weak |28|] . But in the region of phase 
space where the gauge interaction is strong a more refined analysis is needed. It has been 


seen that in this region the effective potential exhibits non-analytic behavior [^]. This issue 
warrants further study. 
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Figure Captions 

Figure 1: Fixed point values of the coupling constants g± and ^2 as N is varied from infinity 

where (gl, g%) = (0, 0) to its lower critical value which is N* = 8 for Np = 4 (bottom curve), 

N* = 4 for Np = 8 (middle curve) and N* = 22 for Np = 80 (top curve). 
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